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ABSTRACT 

Context. There are competing scenarii for planetary systems formation and evolution trying to explain how hot Jupiters came to be 
so close to their parent star. Most planetary parameters evolve with time, making distinction between models hard to do. It is thought 
the obliquity of an orbit with respect to the stellar rotation is more stable than other parameters such as eccentricity. Most planets, to 
date, appear aligned with the stellar rotation axis; the few misaligned planets so far detected are massive ( > 2 Mj). 
Aims. Our goal is to measure the degree of alignment between planetary orbits and stellar spin axes, to detect potential correlation 
with eccentricity or other planetary parameters and to measure long term radial velocity variability indicating the presence of other 
bodies in the system. 

Methods. For transiting planets, the Rossiter-McLaughlin effect allows the measurement of the sky-projected angle f} between the 
stellar rotation axis and a planet's orbital axis. Using the HARPS spectrograph, we observed the Rossiter-McLaughlin effect for six 
transiting hot Jupiters found by the WASP consortium. We combine these with long term radial velocity measurements obtained with 
CORALIE. We used a combined analysis of photometry and radial velocities, fitting models with a Markov Chain Monte Carlo. After 
obtaining ft we attempt to statistically determine the distribution of the real spin-orbit angle if/. 

Results. We found that three of our targets have /? above 90°: WASP-2b: /3 = 153°% WASP-15b: /3 = 139.6°^ 2 and WASP-17b: 
/? = 148.5°^; the other three (WASP-4b, WASP-5b and WASP-18b) have angles compatible with 0°. There is no dependence 
between the misaligned angle and planet mass nor with any other planetary parameter. All orbits are close to circular, with only one 
firm detection of eccentricity on WASP- 18b with e = 0.00848+g{J™j. No long term radial acceleration was detected for any of the 
targets. Combining all previous 20 measurements of /? and our six and transforming them into a distribution of if/ we find that between 
about 45 and 85 % of hot Jupiters have if/ > 30°. 

Conclusions. Most hot Jupiters are misaligned, with a large variety of spin-orbit angles. We find observations and predictions using 
the Kozai mechanism match well. If these observational facts are confirmed in the future, we may then conclude that most hot Jupiters 
are formed from a dynamical and tidal origin without the necessity to use type I or II migration. At present, standard disc migration 
cannot explain the observations without invoking at least another additional process. 

Key words, binaries: eclipsing - planetary systems - stars: individual: WASP-2, WASP-4, WASP-5, WASP-15, WASP-17, WASP-18 
- techniques: spectroscopic 



1. Introduction 

The formation of close-in gas giant planets, the so-called hot 
Jupiters, has been in debate since the discovery of the first of 
them, 51 Pegb, by Mayor & Queloz ( 1995 i. The repeated ob- 



servations of these planets in radial velocity and the discovery 



Send offprint requests to: Amaury. Triaud @unige.ch 

* using observations with the high resolution echelle spectrograph 
HARPS mounted on the ESO 3.6m (under proposals 072.C-0488, 
082.C-0040 & 283.C-5017), and with the high resolution echelle spec- 
trograph CORALIE on the 1.2 m Euler Swiss Telescope, both installed 
at the ESO La Silla Observatory in Chile. The data is made publicly 
available at CDS - Strasbourg 



with HD 209548b (Cha rbonneau et al.|2000| [Henry et al.|2 000) 
that some of them transit has produced a large diversity in plane- 
tary parameters, such as separation, mass, radius (hence density) 
and eccentricity. Although more than 440 extrasolar planets have 
been discovered, of which more than 70 are known to transit, we 
are still increasing the range of parameters that planets occupy; 
diversity keeps growing. 

While it is generally accepted that close-orbiting gas-giant 
planets do not form in-situ, their previous and subsequent evolu- 
tion is still mysterious. Several processes can affect the planet's 
eccentricity and semi-major axis. Inward migration via angular 
momentum exchange with a gas disc, first proposed in|Lin et ah] 
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( |1996| l from work by Goldrei ch & Tremaine| ( |1980] l, is a natural 
and widely-accepted explanation for the existence of these hot 
Jupiters. 

Migration alone does not explain the observed distribu- 
tions in eccentricity and semi-major axis that planets occupy. 
Alternative mechanisms have therefore been proposed such 
as the Kozai mechanism (Kozai 1962; Eggleton & Kiseleva- 



Table 1. List of Observations. The date indicates when the first 
point of the Rossiter-McLaughlin sequence was taken. 



Eggleton 



& Ford 



2001 Wu & Murray 2003]) and planet scattering (jRasio 
996). These mechanisms can also cause a planet to 



migrate inwards, and may therefore have a role to play in the 
formation and evolution of hot Jupiters. These different mod- 
els each predict a distribution in semi-major axis and eccentric- 
ity. Discriminating between various models is done by matching 
the distributions they produce to observations. Unfortunately this 
process does not take into account the evolution with time of the 
distributions and is made hard by the probable combination of a 
variety of effects. 

On transiting planets, a parameter can be measured which 
might prove a better marker of the past history of planets: (5, the 
projection on the sky of the angle between the star's rotation axis 
and the planet's orbital axis. It is believed that the obliquity (the 
real spin-orbit angle if/) of an orbit evolves only slowly and is not 
as much affected by the proximity of the star as the eccentric- 
ity ( ]Hut|1981| |Winn et al.|2005] |Barker & Ogilvie|2009| ). Disc 
migration is expected to leave planets orbiting close to the stel- 
lar equatorial plane. Kozai cycles and planet scattering should 
excite the obliquity of the planet and should provide us with a 
planet population on misaligned orbits with respect to their star's 
rotation. 

As a planet transits a rotating star, it will cause an overall red- 
shifting of the spectrum if it covers the blue-shifted half of the 
star and vice-versa on the other side. This is called the Rossiter- 
McLaughin effect ( |Rossiter| 1 924| |McLaughlin| 1 924\ . It was first 
observed for a pl anet by |Queloz et al.| (|2000[ ). Several papers 
model this e ffect: |Ohta~et al.| < |2005) ; |Gimene z (2006); [GaudT&| 
|WmrHd2007) . 

Among the 70 or so known transiting planets discovered 
since 2000 by the huge effort sustained by ground-based transit- 
ing planet searches, the Rossiter-McLaughlin (RM) effects have 
been measured for 20, starting with observations on HD 209458 
by Queloz et al. ( 2000| ). This method has proven itself reliable at 
giving precise and accurate measurement of the projected spin- 
orbit angle with its best determination done for HD 189733b 
( Triaud et al.| |2009). Basing their analysis on measurements of 
jS in 1 1 systems, 10 of which are coplanar or nearly so, Fabrycky 
( 2009[ ) concluded that the angle distribution is likely 



& Winn 



to be bimodal with a coplanar population and an isotropically 
misaligned population. At that time, the spin-orbit misalignment 
of XO-3b (Hebrard et al. 2008 ) comprised the only evidence of 
the isotropic population. Since then, the misalignment of XO- 
3b has been confirmed by Winn et al. (2009c), and significant 
misalignments have been found for HD 80606b (Moutou et al. 
2009] ) and WASP- 14b ^Johnson et aL]p009] ). Moreover, retro- 
grade orbital motion has been identified in HAT-P-7b (Winn 
et al. 2009b| |Narita et al. 2009). Other systems show indica 



tions of misalignment but need confirmation. One such object is 
WASP- 17b ( |Anderson et al.|2010p which is one of the subjects 
of the present paper. 

The Wide Angle Search for Planets (WASP) project aims at 
finding transiting gas giants (Polla cco etal.|2 006 ). Observing the 
northern and southern hemispheres with sixteen 1 1 cm refractive 
telescopes, the WASP consortium has published more than 20 



Target 


Date 


Instrument 


Paper 


WASP- 18b 


2008/08/21 


HARPS 


this paper 


WASP-8b 


2008/10/05 


HARPS 


|Queloz et al.|{2010| 


WASP-6b 


2008/10/07 


HARPS 


|Uillon et aEfl2009a) 


WASP-4b 


2008/10/08 


HARPS 


this paper 


WASP-5b 


2008/10/10 


HARPS 


this paper 


WASP-2b 


2008/10/15 


HARPS 


this paper 


WASP- 15b 


2009/04/27 


HARPS 


this paper 


WASP- 17b 


2009/05/22 


CORALIE 


this paper 


WASP- 17b 


2009/07/05 


HARPS 


this paper 



transiting planets in a large range of period, mass and radius, 
around stars with apparent magnitudes between 9 and 13. The 
planet candidates observable from the South are confirmed by a 
large radial-velocity follow-up using the CORALIE high reso- 
lution echelle spectrograph, mounted on the 1.2 m Eider Swiss 
Telescope, at La Silla, Chile. As part of our efforts to understand 
the planets that have been discovered, we have initiated a sys- 
tematic program to measure the Rossiter-McLaughlin effect in 
the planets discovered by the WASP survey, in order to measure 
their projected spin-orbit misalignment angles /?. 

In this paper we report the measurement of ft in six south- 
ern transiting planets from the WASP survey, and analyse their 
long term radial velocity behaviour. In sections 2 and 3 we de- 
scribe the observations and the methods employed to extract and 
analyse the data. In section 4 we report in detail on the Rossiter- 
McLaughlin effects observed during transits of the six systems 
observed. In sections 5 and 6 we discuss the correlations and 
trends that emerge from the study and their implications for plan- 
etary migration models. 

2. The Observations 

In order to determine precisely and accurately the angle /?, we 
need to obtain radial velocities during planetary transits at a high 
cadence and high precision. We therefore observed with the high 
resolution echelle spectrograph HARPS, mounted at the La Silla 
3.6 m ESO telescope. The magnitude range within which planets 
are found by the SuperWASP instruments allows us to observe 
each object in adequate conditions. For the main survey pro- 
posal 082.C-0040, we selected as targets the entire population 
of transiting planets known at the time of proposal submission 
to be observable from La Silla during Period 82, i.e. WASP-2b, 
4b, 5b, 6b, 8b and 15b. The result s for WASP-6b are p resented 
separately b y |Gillon et aL| ( |2009a| ) and for WASP-8b by |Queloz 
|et al.| ( [2010[ ). Two targets were added in separate proposals. A 
transit of WASP- 18b was observed during GTO time (072C- 
0488) of the HARPS consortium allocated to this planet because 
of its short and eccentric orbit. During the long-term spectro- 
scopic follow-up of WASP- 17b undertaken for the discovery pa- 
per ( |Anderson et al.|2010| ), three CORALIE measurements fell 
during transit showing a probably retrograde orbit. Observations 
of the Rossiter-McLaughlin with CORALIE confirmed the con- 

pOlO] ) 



elusions of Anderson et al. 



and a follow-up DDT pro- 
posal (283.C-5017) was awardedTime on HARPS. 

The strategy of observations was to take two high precision 
HARPS points the night before transit and the night after transit. 
The radial- velocity curve was sampled densely throughout the 
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Table 2. Stellar parameters used in our model fitting. The v sin / (stellar spectroscopic rotation broadening) and stellar mass estimates 
are used as priors in the analysis. is the micro turbulence. V maC ro is the macrotrubulence. 



Parameters 


units 


WASP-2 (a,b) 


WASP-4 (c) 


WASP-5 (c) 


WASP- 15 (d) 


WASP- 17 (a,e) 


WASP- 18 (f) 


Spectral Type 




Kl 


G8 


G5 


F7 


F4 


F6 


Teff 


K 


5150 ± 80 


5500 ± 100 


5700 ± 100 


6300 ± 100 


6650 ± 80 


6400 ± 100 


B—V 


mag 


0.86 ± 0.05 


0.73 ± 0.05 


0.66 ± 0.05 


0.48 ± 0.05 


0.38 ± 0.05 


0.45 ± 0.05 


logg 


dex 


4.40 ±0.15 


4.5 ± 0.2 


4.5 ± 0.2 


4.35 ±0.15 


4.45 ±0.15 


4.4 ±0.15 


[Fe/H] 


dex 


-0.08 ± 0.08 


-0.03 ± 0.09 


+0.09 ± 0.09 


-0.17 ± 0.11 


-0.19 ± 0.09 


0.00 ± 0.09 


logR' HK 


dex 


-4.84 + 0.10 


-4.50 ± 0.06 


-4.72 ± 0.07 


-4.86 ± 0.05 




-4.85 ± 0.02 


ft 


kms 1 


0.9 ± 0.1 


1.1 ± 0.2 


1.2 ±0.2 


1.4 ±0.1 


1.7 ± 0.1 


1.6 ±0.1 


^macro 


kms 1 


1.6 ±0.3 


2.0 ± 0.3 


2.0 ± 0.3 


4.8 ± 0.3 


6.2 ± 0.3 


4.8 ± 0.3 


v sin / 


kms" 1 


1.6 ±0.7 


2.0 ± 1.0 


3.5 ± 1.0 


4.0 ± 2.0 


9.8 ± 0.5 


11.0 ± 1.5 




M 


0.84 ±0.11 


0.93 ± 0.05 


1.00 ±0.06 


1.18 ± 0.12 


1.2 ±0.12 


1.24 ±0.04 



Notel. refe rences: (a) this paper, (b) |Cameron etlLlpOOT) , (c) |Gillon et aL] ( |2009c^ (d) |West et al"7|p009) , (e) |Anderson"etaT1p0T0) ; (f) |Hellier| 
[eTall ( |2009l 



transit, beginning 90 minutes before ingress and ending 90 min- 
utes after egress. The data taken before ingress and after egress 
allow any activity-related offset in the effective velocity of the 
system's centre of mass to be determined for the night of obser- 
vation. In addition, radial velocity data from the high resolution 
echelle spectrograph CORALIE mounted on the Swiss 1.2 m 
Euler Telescope, also at La Silla was acquired to help search for 
a long term variability in the the periodic radial velocity signal. 

All our HARPS observations have been conducted in the 
OBJO mode, without simultaneous Thorium-Argon spectrum. 
The velocities are estimated by a Thorium-Argon calibration at 
the start of the night. HARPS is stable within lms" 1 across a 
night. This is lower than our individual error bars and leads to 
no contamination of the Th-Ar lamp onto the stellar spectrum 
easing spectral analysis. 



3. The Data Analysis 

3.1. Radial-velocity extraction 

The spectroscopic data were reduced using the online Data 
Reduction Software (DRS) which comes with HARPS. The ra- 
dial velocity information was obtained by removing the instru- 
mental blaze function and cross-correlating each spectrum with 
one of two masks. This correlation is compared with the Th-Ar 
spectrum acting as a reference; see |Baranne et al. ( 1996}, Pepe 
etaLlp002) & |Mayor et al^pPOl] ) for details Recently the DRS 
was shown to achieve remarkable precision (Mayor et al. 2009) 
thanks to a revision of the reference lines for Thorium and Argon 
by |Lovis & Pep e (2007 ). Stars with spectral type earlier than G9 
were reduced using the G2 mask, while those of K0 or later were 
cross-correlated with the K5 mask. A similar software package 
is used for CORALIE data. A resolving power R = 110000 
for HARPS yields a cross-correlation function (CCF) binned in 
0.25 km s 1 increments, while for CORALIE, with a lower res- 
olution of 50000, we used 0.5 km s" 1 . The CCF window was 
adapted to be three times the size of the full width at half maxi- 
mum (FWHM) of the CCF. 

All our past and current CORALIE data on the stars pre- 
sented here were reprocessed after removal of the instrumen- 
tal blaze response, thereby changing slightly some radial veloc- 



ity values compared to those already published in the literature. 
Correcting this blaze is important for extracting the correct RVs 
for the RM effect. The uncorrected blaze created a slight sys- 
tematic asymmetry in the CCF that was translated into a bias in 
radial velocities. 

1 cr error bars on individual data points were estimated from 
photon noise alone. HARPS is stable long term within lms" 1 
and CORALIE at less than 5ms" 1 . These are smaller than our 
individual error bars and thus have not been taken into account. 



3.2. Spectral analysis 

Spectral analysis is needed to determine the stellar atmospheric 
parameters from which limb darkening coefficients can be in- 
ferred. We carried out new analyses for two of the target stars, 
WASP-2 and WASP-17, whose previously-published spectro- 
scopic parameters were of low precision. For our other targets, 
the atmospheric parameters were taken from the literature, no- 
tably the stellar spectroscopic rotation broadening v sin /[[] 

The individual HARPS spectra can be co-added to form an 
overall spectrum above S/N ~ 1 : 100, suitable for photospheric 
analysis which was performed using the uclsyn spectral synthe- 
sis package ( Smith|1992 Smalley et al.|200T i and atlas9 mod- 
els without convective overshooting (Castelli et al. 1997 1 and the 
same method as described in many discovery papers published 
by the WASP consortium (eg: [Wilson et al.|(|2008)). 



The stellar rotational v sin / is determined by fitting the pro- 
files of several unblended Fe i lines. The instrumental FWHM 
was determined to be 0.065A from the telluric lines around 
6300A. 

For WASP-2, a value for macrotu rbulence (v mac ) of 1.6 ± 
0.3 kms -1 was adopted (Gray 20081. A best fitting value of 
vsin/ = 1.6 + OJkms" 1 was obtained. On WASP-17, a value 
for macroturbulence (v mac ) of 6.2 + 0.3 kms" 1 was used (Gray 
2008). The analysis gives a best fitting value of vsin/ = 9.8 + 



1 throughout this paper we use the symbol / to denote the inclination 
of the stellar rotation axis to the line of sight, while i represents the 
inclination of the planet's orbital angular momentum vector to the line 
of sight 
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0.5 km s 1 . The error on v mac is taken from the scatter around fit 
Gray (2008 ) and is propagated to the v sin / . 



to 



B—V were estimated from the effective temperature and used 
in the calculations of the logR' HK dNoyes et al. 1984 Santos et al. 
2000 Boisse et aH2 009). Errors refer to the photon noise: they 



do not include systematic effects likely to arise and affect low 
values of log R' HK due to the low signal to noise in the blue or- 
ders. WASP-17 does not have a value since this stellar activity 
indicator is only calibrated for B - V e [0.44, 1.20]. 

All stellar parameters, used as well as derived, are presented 
in Table |2] 

3.3. Model fitting 

The extracted radial velocity data was fitted simultaneously with 
the transit photometry available at the time of analysis. Three 
models are adjusted to the data: a Keplerian radial velocity or- 
bitjHilditch 2001), a photometric planetary transit ( |Mandel &| 
|Agol|2002| l, and a spectroscopic transit, also known as Rossiter- 
McLaughlin effect (Gimenez 2006). This combined approach is 
very useful for taking into account all of the possible contribu- 
tions to the uncertainties due to correlations among all relevant 
parameters. A single set of parameters describes both the pho- 
tometry and the radial velocities. We use a Markov Chain Monte 
Carlo (MCMC) approach to optimize the models and estimate 
the uncertainties of the fitted parameters. The fit of the model to 
the data is quantified using the^ 2 statistic. 

The code is described in detail by |Triaud et al. (20091), has 
been used several times (eg:|Gillon et al.|(|2009a )) and is similar 
to the code described in Cam eron et al.|P007| ). 

We fitted up to 10 parameters, namely the depth of the pri- 
mary transit D, the radial velocity (RV) semi-amplitude K, the 
impact parameter b, the transit width W, the period P, the epoch 
of mid-transit To, e cos oj, e sin a>, V sin / cos B, and V sin / sin B. 
Here e is the eccentricity and u the angle between the line of 
sight and the periastron, V sin / is the sky-projected rotation ve- 
locity of the statj^jwhile B is the sky-projected angle between the 
stellar rotation axis (Hosokawa 1 953[ [Giiri enez 2006) and the 
planet's orbital axi^] 

These parameters have been chosen to reduce correlations 
between then. The use of uncorrected parameters allows to 
explore parameter space more efficiently since the correlation 
length between jumps is smaller. Eccentricity and periastron an- 
gle were paired as were V sin / and B. This breaks a correlation 
between them (the reader is invited to compare Figs.|2jl &[3]for 
a clear illustration for choosing certain jump parameters as op- 
posed to others). This way we also explore solutions around zero 
more easily: e cos a> and e sin a> move in the ] - 1 , 1 [ range while e 
could only be floating in ]0,1 [. For exploring particular solutions 
such as a circular orbit, parameters can be fixed to certain values. 

We caution that, as noted by |Ford| (12006 ) that the choice 
of ecosw and esinw as jump variables implicitly imposes a 
prior that is proportional to e. This approach thus has a ten- 
dency to yield a higher eccentricity than would be obtained with 
a uniform prior, in cases when e is poorly-constrained by the 
data. A similar argument applies to the use of V sin / cos B and 



2 we make a distinction between v sin / and V sin /: v sin / is the value 
extracted from the spectral analysis, the stellar spectroscopic rotation 
broadening, while V sin / denotes the result of a Rossiter-McLaughlin 
effect fit. Both can at times be different. Each, although caused by the 
same effect, is independently measured making the distinction worth- 
while. 

3 y8 = -A, another notation used in the literature for the same angle. 



V sin / sin B as jump variables, in cases where the impact param- 
eter is low and there is a strong degeneracy between V sin / and 
B in modelling the Rossiter-McLaughlin effect. In such cases, 
however, the tendency to overestimate V sin / from the Rossiter- 
McLaughlin effect can effectively be curbed by imposing an in- 
dependent, spectroscopically-determined v sin / prior on V sin /, 
as we have done here. 

In addition to the physical free floating parameters, we need 
to use one y velocity for each RV set and one normalisation 
factor for each lightcurve as adjustment parameters. These are 
found by using optimal averaging and optimal scaling, y veloc- 
ities represent the mean radial velocity of the star in space with 
respect to the barycentre of the Solar System. Since our analysis 
had many datasets, the results for these adjustment parameters 
have been omitted, not adding anything to the discussion. 

During these initial analyses we also fitted an additional ac- 
celeration in the form of an RV drift y but on no occasion was 
it significantly different from zero. We therefore assumed there 
was no drift for any of our objects. We will give upper limits for 
each star in the following sections. 

The MCMC algorithm perturbs the fitting parameters at each 
step i with a simple formula: 



P i j=P i -i,j+fo- Pj G(0,l) 



(1) 



where Pj is a free parameter, G is a Gaussian random number of 
unit standard deviation and zero mean (meaning a Gaussian prior 
on each parameter), while <x is the step size for each parameter. 
A factor / is used to control the chain and ensures that 25 % of 
steps are being accepted via a Metropolis-Hastings algorithm, 



as recommended in |Tegmark et al. ( 2004 ) to give an optimal 



exploration of parameter space. 

The step size is adapted by doing several initial analyses. 
They are adjusted to produce as small a correlation length as 
possible. Once the value is chosen, it remains fixed. Only / fluc- 
tuates. 

A burn-in phase of 50000 accepted steps is used to make 
the chain converge. This is detected when the correlation length 
of each parameter is small and that the average x 1 does not im- 
prove anymore (Tegmark et al. 2004). Then starts the real chain, 
of 500 000 accepted steps, from which results will be extracted. 
This number of steps is used as a compromise between computa- 
tion time and exploration. Statistical tests, notably by comparing 
X 2 are used to estimate the significance of the results. 

Bayesian penalties acting as prior probability distribution 
can be added to x 1 to account for any prior information that 
we might have on any fitted or derived parameter. Stellar mass 
M+ can notably be inserted via a prior in the MCMC in order to 
propagate its error bars on the planet's mass. We also inserted the 
v sin / found by spectral analysis as priors in some of our fits to 
control how much the fit was dependent on them and the result- 
ing value of V sin / and whether this influenced the fitted value 
of B. Because of the random nature of an MCMC, sometimes a 
step with an impact parameter close to zero is taken. This can 
cause V sin / to wander to unphysical values because of the de- 
generacy between V sin / and B at low impact parameters. This 
is controlled by imposing a prior. The prior values are in Table 

El 

We use a quadratic limb-darkening law with fixed values 
for the two limb darkening coefficients appropriate to the stel- 
lar effective temperature. They were extracted for the photom- 
etry from tables published in Clare t| ([2000 ). For the radial ve- 
locity (the Rossiter-McLaughlin effect is also dependent on limb 
darkening) we use values for the V band. Triaud et al. (2009 ) 
showed that HARPS is centred on the V band. The coefficients 
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Fig. 1. Fit results for WASP-2b. a) Overall Doppler shift reflex motion of the star due to the planet and residuals, b) Zoom on the 
Rossiter-McLaughlin effect and residuals. Black inverted triangles are SOPHIE data, black triangles represent CORALIE points, 
red dots show the HARPS data. The best fit model is also pictured as a plain blue line. In addition to our best model found 
with Vsinl = 0.99 km s~' we also present models with no RM effect plotted as a dotted blue line, RM effect with /3 — and 

V sin / = 0.9 km s _1 drawn with a dashed-dotted blue line and RM effect with /3 — and V sin / = v sin I =1.6 km s pictured with 
a dashed-double dotted blue line. In the residuals, the open symbols represent in the values with the size of the circle decreasing 
with the likelyhood of the model, c) Posterior probability distribution issued from the MCMC showing the distribution of points 
between e cos a> and e sin m.d) Posterior probability distribution issued from the MCMC showing the distribution of points between 

V sin / cos/? and V sin / sin/3. The black disc shows where the distribution would be centred only changing to ft = 0. The dotted 
line shows where zero is. The straight lines represent the median of the distribution, the dashed lines plot the position of the average 
values, the dash-dotted lines indicate the values with the lowest^ 2 (some lines can overlap). The size of boxes c) and d) represents 
7 times the 1 <x distance on either side of the median. 



were chosen for atmospheric parameters close to those presented 
in Table |2] 



3.4. Extracting the results 

For each star, we performed four analyses, each using a MCMC 
chain with 500 000 accepted steps: 

- 1 . a prior is imposed on V sin /, eccentricity is fixed to zero; 

- 2. no prior on V sin /, eccentricity is fixed to zero; 

- 3. a prior is imposed on V sin /, eccentricity is let free; 

- 4. no prior on V sin/, eccentricity is let free. 

This is to assess the sensitivity of the model parameters to a 
small but uncertain orbital eccentricity and to the v sin / value 



found by spectral analysis which, as demonstrated in Triaud 
|et al.| ( |2009| l, can seriously affect the fitting of the Rossiter- 
McLaughlin effect. The comparative tables holding the results of 
these various fits are available in the appendices to support the 
conclusions we reach while allowing readers to form their own 
opinion. In addition, we also conducted control chains fixing 
the parameters controlled by the photometry in order to check 
whether this was a limiting source of errors in the determina- 
tion of our most important parameter: ft. The results from these 
chains are in the appendices as well. Although different in their 
starting hypotheses all the chains are also useful at checking 
their respective convergence. Our final results are presented in 
Table [3] 
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The best solution is found in the best of the four fits by com- 
paring^- 2 and using Ockham's principle of minimising the num- 
ber of parameters for similar results: for fits with similar v 2 . , 

1 ™ reduced 

we usually choose a circular solution with no prior on V sin /. 
Results are extracted from the best fit by taking the median of 
the posterior probability distribution for each parameter, deter- 
mined from the Markov chain. Errors bars are estimated from 
looking at the extremes of the distribution comprising the 68.3 % 
confidence region of the accepted steps. The best solution is not 
taken from the lowest^ 2 as it is dependent on the sampling and 
chance encounter of a - small - local minimum. Scatter plots will 
be presented with the positions of the best x 2 , the average and 
the median for illustration. 

In the following section and in tables, several statistical val- 
ues are used: x 2 is the value found for all the data, while x\s/ 
gives the value ofx 2 solely for the radial velocities. The reduced 
X 1 for the radial velocities, denoted by ;^ educed , is used to esti- 
mate how well a model fits the data and to compare various fits 
and their respective significance. In addition we will also use the 
residuals, denoted as O - C. These estimates are only for radial 
velocities. The results from photometry are not mentioned since 
they are not new. They are only here to constrain the shape of 
the Rossiter-McLaughlin effect. 

When giving bounds, for eccentricity and long term radial 
velocity drift, we quote the 95% confidence interval for exclu- 
sion. 



4. The Survey Results 

4.1. WASP-2b 

A sequence of 26 RV measurements was taken on WASP-2 us- 
ing HARPS on 2008 October 15, with additional observations 
made outside transit as given in the journal of observations pre- 
sented in the appendices. The cadence during transit was close 
to a point every 430s. The average photon noise error of that 
sequence is 5.7 ms -1 . We made additional observations with 



CORALIE to refine the orbital solution obtained by Cameron 
et al.| ( |2007) using the SOPHIE instrument on the 1.93 m tele- 
scope at Observatoire de Haute-Provence, and to look for long- 
term variability of the orbit. 20 measurements were taken with 
a mean precision of 13.9 m s _1 over close to 1 1 months between 
2008 October 25 and 2009 September 23. All the RV data is 
available in the appendices along with exposure times. 

To establish the photometric ephemeris and the transit geom- 
etry, we fitted the photometric datasets of Cameron et al. ( 2007| > 
(3 seasons by SuperWASP in the unfiltered WASP bandpass), 
Charbonneau et al. (2007 1 (a z band Keplercam lightcurve) and 
Hrudkova et al. (2009) (a William Herschel Telescope AG2 R 
band transit curve). 

WASP-2b's data were fitted with up to 10 free parameters 
plus 8 independent adjustment parameters: three y velocities for 
the three RV data sets and five normalisation factors for photom- 
etry. This sums up to 58 RV measurements and 8951 photometric 
observations. 

^reduced does not i m P rove significantly between circular and 
eccentric models. We therefore impose a circular solution. The 
presence of a prior on V sin / does not affect the results. We find 
Vsin/ = 0.99 



-o 



found in section 



" 7 kms 1 in accordance with the vsin /value 
3.2] The fit delivers p = 153°+^. The overall 
root-mean-square (KMS) scatter of the spectroscopic residuals 
about the fitted model is 1 1 .73 m s -1 . During the HARPS transit 
sequence these residuals are at 6.71 ms _1 . 



Fig. [3] shows the resulting distribution as Vsin I vs. p. We 
detect V sin / significantly above zero with confidence interval 
showing that 99.73% (3 cr) of the posterior probability function 
has Vsin / > 0.2kms _1 while p > 77.26°. We have computed 6 
additional chains in order to test the strength of our conclusions. 



Table A.l shows the comparison between the various fits; we 



invite the reader to refer to it as only important results are given 
in the text. 

In all cases, eccentricity is not detected being below a 3 cr 
significance from circular which is likely affected from the poor 
coverage of the phase by the HARPS points. Circular solutions 
are therefore adopted. We fix the eccentricity's upper limit to 
e < 0.070. In addition no significant long term drift was detected 
in the spectroscopy: \y\ < 36 ms -1 yr _1 . 

Using the spectroscopically-determined vsin/ value of 
1.6kms _1 and forcing p to zero, AfLjuced changes from2.14±0.27 
to 3.49 ± 0.39, clearly degrading the solution. We are in fact 
7.6 cr away from the best-fitting solution, therefore excluding an 
aligned system with this large a V sin /. This is also excluded by 
comparison to a fit with a flat RM effect at the 6.7 cr. Similarly, 
a fit with an imposed Vsin/ = 0.9 kms 1 and aligned orbit is 
found 5.6 cr from our solution. On Fig[lJ), we have plotted the 
various models tested and their residuals so as to give a visual 
demonstration of the degradation for each of the alternative so- 
lutions. 



4.2. WASP-4b 

We obtained a RM sequence of WASP-4b with HARPS on 2008 
October 8; other, out of transit, measurements are reported in 
the journal of observations given in the appendices. The RM 
sequence comprises 30 data points, 13 of which are in transit, 
taken at a cadence of 630 s with a mean precision of 6.4 m s . 
The spectrograph CORALIE continued monitoring WASP-4 and 
we add ten radial velocity measurements to the ones published 
in | Wilson et al.| ( |2008| ). These new data were observed around 
the time of the HARPS observations, about a year after spectro- 
scopic follow-up started. 

In photometry we gathered 2 timeseries in the WASP band- 
pass from Wil son et al.| ( |2008[ ) and an R band C2 Euler transit 
plus a VLT/FORS2 z band lightcurve obtained from Gillon etaT] 
(2009c ) to establish the transit shape and timing. 

The WASP-4b data were fitted with up to 10 free parameters 
to which 6 adjustment parameters were added: two y velocities 
for RVs and four normalisation factors for the photometry. In to- 
tal, this represents 56 radial velocity points and 9989 photomet- 
ric measurements. Gillon et al. (2009c) let combinations of limb 
darkening coefficients free to fit the high precision VLT curve. 
We used and fixed our coefficients on their values. 

Because the impact parameter is small, a degeneracy be- 
tween /3 and Vsin/ appeared, as expected (see Figs. [2}i &|3). 
The values on stellar rotation for our unconstrained fits reach 
unphysical values as high as Vsin/ = 150kms _1 . We imposed 
a prior on the stellar rotation to restrict it to values consistent 
with the spectroscopic analysis. 

The reduced x 1 is the same within error bars whether ec- 
centricity if fitted or fixed to zero. Therefore the current best 
solution, by minimising the number of parameters, is a circular 
orbit. 

The eccentricity is constrained to e < 0.0182. Thanks to the 
long time series in spectroscopy we also investigated the pres- 
ence of a long term radial velocity trend. Nothing was signifi- 



cantly detected: \y\ < 30 m s yr 
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Fig. 2. Fit results for WASP-4b. Nota Bene: Legend similar to the legend in PigJT] 



Because of the small impact parameter the spin-orbit angle 
is poorly constrained with B = -4°^, even when a prior is 
imposed on Vsin/. The high S/N of the Rossiter-McLaughlin 
effect allows us to exclude a projected retrograde orbit. 



4.3. WASPSb 

Using HARPS, we took a series of 28 exposures on WASP-5 at 
a cadence of roughly 630s with a mean photon noise of 5.5 m s" 1 
on 2008 October 16. Other measurements were obtained at dates 
before and after this transit. Five additional CORALIE spectra 
were acquired the month before the HARPS observations. They 
were taken about a year after the data published in |Anderson 



et al. ( 2008 1. All spectroscopic data is available from the appen- 



dices. 

To help determine transit parameters, published photometry 
was assembled and comprises three seasons of WASP data, two 
C2 Euler lightcurves in R band, and one FTS i' band lightcurve 
( |Anderson et al.|2008| >. 

WASP-5b's 49 RV measurements and 14754 photometric 
points were fitted with up to 10 free parameters to which 8 ad- 
justment parameters had to be added: two y velocities and six 
normalisation factors. 



The imposition of a prior on Vsin/ prior makes little im- 
pact on the final results (see appendices) and their lcr error bars 
but prevents Vsin/ to go to unphysical values when, through 
the random process of the MCMC, the impact parameter b gets 
very close to zero on a few occasions. We choose the solution 
using a prior. The priorless solution gives a V sin / fully consis- 
tent with v sin / thereby obtaining an independent measurement 
of the projected stellar equatorial rotation speed. Allowing ec- 
centricity to float did not produce a significantly better fit. It has 
a 99.6 % chance of being different from zero: at 2.9 <x. Thus, 
minimising the number of parameters for a similar fit, we chose 
the solution with a circular orbit and simply place an upper limit 
on the eccentricity: e < 0.0351. No long term RV trend appears 
at this date: \j\ < 47 m s~' yr" 1 . 

Parameters extracted are similar to those that were pub- 
lished [jnjGnion[eFlI]d200^ & |Anderson eFaT] pOOg) and 
with Southworth et al. ( 2009 1 using a independent dataset. 



The projection of the spin-orbit angle is found to be: B = 
-12.1°^g g and we obtain an independent measurement of 



Vsin/ — 3.24^35 km s fully compatible with the spectral 
value that was used as a prior in other fits. Results are presented 
in Table [3] 

The^ , i 2 educed for spectroscopy (see Table A. 1 1 is quite large, at 
3.68 ±0.44. The O-C for CORALIE data stand at 17.94 ms" 1 to 
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Fig. 4. Fit results for WASP-5b. Nota Bene: Legend similar to the legend in FigJT] 



be compared with an average error bar of 18.13 ms -1 . The bad- 
ness of fit therefore comes from the HARPS sequence which has 
a dispersion of 8.98 m s -1 for an average error bar of 5.49 ms -1 . 
From Fig.|4|) we can see that residuals are quite important during 
the transit; Fig.|4fl also shows that the MCMC does not produce 
a clean posterior distribution. This is mostly caused by impact 
parameter values nearing zero during parameter exploration and 
causing a degeneracy between V sin / & B. This can be observed 
on Figure[5]with similitude to what occurs to WASP-4. 



No better solution can be adjusted to the data: we remind 
that the RM effect is fitted in combination with six photomet- 
ric sets which strongly constrain the impact parameter, depth 
and width of the Rossiter-McLaughlin effect. The V sin / cos B 
vs V sin / sin/3 distribution is not centred on zero but close to it. 
This may come from the intrinsic dispersion in the data. Among 
the six data points which are spread over the rest of the phase, 
we have a dispersion of 11.92ms -1 . A likely cause to explain 
the data dispersion is stellar activity. Table [2] indicates that this 
star is moderately active. A longer discussion on ^ 1 2 educed > 1. is 
presented in section |5.3| |Santos et al.| ( |2000) > show that for the 
log R' HK that we find, we can expect a variation in velocities of 
the order of 7 to 12 m s" 1 . 



4.4. WASP-15b 

Observations were conducted using the spectrographs 
CORALIE and HARPS. 23 new spectra have been ac- 
quired with CORALIE in addition to the 21 presented in [West 
|et al.| ( |2009[ ) and extending the time series from about a year 
to 500 days. We observed a transit with HARPS on 2009 April 
27. 46 spectra were obtained that night, 32 of which are during 
transits with a cadence of 430s. Additional observations have 
been taken as noted in the journal of observations. 

The photometric sample used for fitting the transit has data 
from five time-series in the WASP bandpass, as well as one / and 
one R band transit from C2 Euler (West et al. 2009 ). The spectral 
data were partitioned into two sets: CORALIE and HARPS. 

7 normalisation factors and 2 y velocities were added to ten 
free floating parameters to adjust our models to the data which 
included a total of 95 spectroscopic observations and 23 089 
photometric measurements. 

For the various solutions attempted, v 2 , , are found the 

, . 1 ' ^ reduced 

same (Table A. 2 1. We therefore choose the priorless, circular ad- 



justment as our solution. 

Compared to [West et al.| ([2009), parameters have only 
changed little. Thanks to the higher number of points we give 
an upper limit on eccentricity: e < 0.087 (Fig. [5}; shows re- 
sults consistent with zero); there is no evident long term evolu- 
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tion in the radial velocities, which is constrained within: \y\ < Table 4. List of y velocities for WASP-17's RV sets. 

llms 'yr 1 . The projected spin-orbit angle is found rather 

large with B = 139.6° ^2 making WASP- 15b appear as a retro- 
grade planet with a very clear detection. V sin / is found within 



Instrument Dataset 



y (m s ') 



1 o~ o f the spectrally analysed value of v sin / from West et al. 
(2009) at 4.27+°-^ kms 
dependent measurement 



and as such constitutes a precise in- 



X. 



reduced 



1.51+0.19 for the spectroscopy, indicating a good 
fit of the Keplerian as well as of the Rossiter-McLaughlin effect, 
the best fit in this paper. Full results can be seen in Table [3] 



4.5. WASP-17b 

On 2009 May 22, 1 1 CORALIE spectra were obtained at a ca- 
dence of 2030s with an average precision of 33.67 m s _1 to con- 
firm the detection of retrograde orbital motion announced by 
|Anderson et al.| ( |2010) . The sequence was stopped when airmass 
reached 2. HARPS was subsequently used and on 2009 July 5 
a sequence of 42 spectra was acquired with a cadence of 630s 
during transit. They have a mean precision of 19.02 m s _1 . In ad- 
dition to these and to data already published 12 CORALIE spec- 
tra and 15 HARPS spectra were obtained. All the spectroscopic 
data is presented in the appendices. 



CORALIE Rossiter-McLaughlin effect -49500.80!: 

CORALIE orbital Doppler shift -49513.67!° 

HARPS Rossiter-McLaughlin effect -49490.59! 

HARPS orbital Doppler shift -49491.68 



Ail 



The photometry includes five timeseries of data in the WASP 
bandpass, and one C2 Euler I band transit (Anderson et al. 

|20"T0| 



The model had to adjust up to 10 free floating parameters and 
10 adjustment parameters (6 photometric normalisation factors 
and 4 radial velocity offsets) to 15 690 photometric data points 
and 124 spectroscopic points. 

The RV was separated into four datasets fitted separately as 
detailed in Table [4] This was done to mitigate the possibility 
that the RM effect was observed at a particular activity level for 
the star. Stellar activity adds an additional RV variation. For a 
set where this data is taken randomly over some time, one ex- 
pects activity to act like a random scatter around a mean which 
would be the true y velocity of the star in space. But for a se- 
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Fig. 6. Fit results for WASP-17b. On a) and b) black circles represent the RM effect taken with CORALIE, while black triangles 
picture the remaining CORALIE measurements; red dots show the HARPS RM data, red triangles are the remained HARPS points. 
Nota Bene: Legend similar to the legend in Fig{T] 



quence such as the RM effect, we expect only a slowly-varying 
radial-velocity bias caused by the activity level on the star on 



the night concerned. This analysis method is explained in Triaud 
|et al.| ( |2009| l which showed an offset in y velocities between 
different Rossiter-McLaughlin sequences of HD 189733 which 
can only be attributed to stellar variability. The large number of 
CORALIE and HARPS measurements outside transit and their 
large temporal span allowed us to separate RV sets for WASP- 17 
but not for the other targets. Table |4] shows the four values of y. 
We remark a difference of 1 3 m s~' for CORALIE, justifying our 
segmentation of the data. 

Among the four computed chains, we select the circular so- 
lution, with prior on V sin / since our results show eccentricity 
is not significantly detected but that the prior on V sin / prevents 
the MCMC from wandering to small impact parameters leading 
to the degeneracy between V sin / and B. 



The non significant eccentricity presented by Anderson et al. 
(2010) was not confirmed, so a circular orbit was adopted. We 
confine to within e < 0.110. Eccentricity affects the derived 
value of the stellar density, and thereby also affects the planet's 
radius measurement. Our circular solution suggests that WASP- 
17b's radius is 1.986^S - [5? ^j, making it the largest and least 



dense extrasolar planet discovered so far. We looked for an addi- 
tional long term acceleration but found none: \y\ < 18 m s~' yr~' . 

The Rossiter-McLaughlin effect is well fitted. The residu- 
als show some dispersion about the model during the HARPS 
sequence. At the end of the HARPS transit, the airmass at- 
tained high values which account for the larger error bars, the 
sparser sampling and higher dispersion. By comparison the 
CORALIE sequence appears better: its longer exposures blurred 
out short-term variability. Both Vsin/ and B are unambigu- 
ously detected. WASP- 17b is on a severely misaligned orbit: 
Vsin/ = 9.9210ns- 1 and B = 148.5°+^. Full results are dis- 
played in Table [3] 



4.6. WASP-18b 

Soon after WASP- 18b was confirmed by the spectrograph 
CORALIE, a Rossiter-McLaughlin effect was observed with 
HARPS. We obtained 19 measurements at a cadence of 630s 
on 2008 August 21. The mean photon noise for the transit se- 
quence is 6.99 ms~'. Seeing and airmass improved during the 
sequence, increasing the S/N and decreasing the individual error 
bars. Additional data were also acquired out of transit. |Hellier| 
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etaE1 ( |2009[ > presented 9 RV measurements from CORALIE. 28 
more have been taken and are presented in this paper. They span 
over three months. The total data timeseries spans close to 500 
days. All RV measurements are presented in the journal of ob- 
servations at the end of the paper. 

Transit timing and geometry were secured by four photomet- 
ric series: two SuperWASP seasons and two C2 Euler transits in 
R band, presented in Hellier et al. ( 2009| l. 



The fitted data comprises 8593 photometric measurements 
and 60 radial velocities. Ten free parameters were used, with, in 
addition, four normalisation constants and two y velocities. 

Eccentricity is clearly detected, improving deduced from 
5.58 ± 0.47 to 3.70 ± 0.36 (from 4.31 ± 0.46 to 2.00 + 0.32 if 
we remove the RM effect from the calculation). We therefore 
exclude a circular solution. 

The V sin / found in the priorless chain differs from the spec- 
tral analysis (15.57^ ^ instead of 1 1 ± 1.5 km s -1 ), this solution 
is preferred so as to not produce biased results. For this par- 
ticular case, we should consider V sin / more like a amplitude 
parameter in order to fit the Rossiter-McLaughlin effect rather 
than a bone fide measurement of projected rotation of the star. 
Therefore, the solution we favour is that of an eccentric orbit, 
without a prior on the V sin /. 

Results are presented in Table. [3] and the best fit is shown in 
Fig. [7] This Rossiter-McLaughlin effect is one of the largest so 



far measured, with an amplitude of nearly 185 m s 1 . During the 
transit sequence O — C — 15.02 ms~' for a mean precision of 

6.95 ms" 1 : the fit is poor; deduced = 3 - 70 ± °- 36 - This is likel y 
caused by a misfit of a symmetric Gaussian on a no longer sym- 
metrical CCfQ We are in fact resolving the planet transit in front 
of the star like spots can be detected via Doppler tomography. 
This has recently been modelled and detected for HD 189733b, 
as a Doppler shadow by |Cameron et al.| ( 2010 >. The accuracy on 
the jS parameter is not affected by the misfit since it is measured 
from the asymmetry of the Rossiter-McLaughlin effect. It is in 
essence estimated from the difference of time spent between the 
two hemispheres of the star. 

Therefore all parameters can be trusted except the Vsin/, 
including the much sought after f3 angle. We find it to be consis- 
tent with zero within 1.5 cr: y3 = -4.0° ^H. The precision on this 
angle is the best we measured, something that is not reflected in 
the fit, we therefore doubled error bars to (3 — -4.0° This 
is in part thanks to the brightness of the star, allowing precise 
measurements of a large amplitude effect. Any departure from 
the model is quickly penalised in x 2 by the data. Similarly, ec- 
centricity is detected above 9 cr with e = 0.00848^ ^"^ thanks 



this was noted in 



Triaud et al. 



2009j in the case o f HP 189733b 



and CoRoT-3b, but can also be seen on fits of CoRoT-2b (Bouch y"et al.| 
[2008) , Hat-P-2b JLoeillet et al.|2008j > and others. 
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to the large amplitude of the reflex motion. Note, that fitting 
e cos U) and e sin u> can correspond to fitting e proportional to e 
and tending to bias the search for solutions towards higher val- 
ues. We attempted a few control fits exploring instead y[e cos u 
& y[e sin lo. The results showed there is no bias in our analy- 
sis, so strongly is the eccentricity constrained by the radial ve- 
locities. The spectroscopic coverage gives us the chance to put 
some limits on an undetected long term radial velocity drift: 
\y\ < 43 ms -1 yr -1 . 

The other parameters are consistent with what has been pub- 
lished by Hellier et al. 2009 and are presented in Table [3] 



5. Overall results 

Our fits to the Rossiter-McLaughlin effect confirm the presence 
of planetary spectroscopic transit signatures in all six systems. 
While three of the six appear closely aligned, the other three 
exhibit highly-inclined, apparently retrograde orbits. The orbits 
of all six appear close to circular. Only the massive WASP-18b 
yields a significant detection of orbital eccentricity. 



5.1. Orbital eccentricities 



As observed in Gillon et al. (2009b i, treating eccentricity as a 
free fitting parameter increases the error bars on other parame- 
ters; we are exploring a larger parameter space. One might argue 
that allowing eccentricity to float is necessary since no orbit is 
perfectly circular, therefore making an eccentric orbit the sim- 
plest model available. We argue against this for the simple rea- 
son that if statistically we cannot make a difference between an 
eccentric and a circular model then it shows that the eccentric 
model is not detected. Actually, the mere fact of letting eccen- 
tricity float biases the result towards a small non zero number, 



a bias which can be larger than the actual physical value (Lucy 
& Sweeney|1971[ l. Hence letting eccentricity float when it is not 
detected is to allow values of parameter space for all parameters 
to be explored which do not need to be. This is why, unless^ 2 is 
significantly improved by adding two additional parameters to a 
circular model, we consider the former as preferable. In addition 
to the risk of biasing, there is a strong assumption that due to 
tidal effects circularising the planet's orbit, eccentricities are re- 
ally small and therefore undetectable for the majority of targets. 
It is therefore reasonable to assume a value of zero when the data 
does not contradict it. To facilitate comparison, we also present 
the results of fits with floating eccentricity. These are given in 
the appendices; our preferred solutions are described in the text 
and in table [3] 

Only for WASP- 18b, have we detected some eccentricity in 
the orbit, thanks primarily to the high amplitude of the RV signal 
and the brightness of the target. The amount of RV data taken 
on WASP-18b is not really more than for the other targets. In 
addition to a high semi-amplitude, sampling is another key to 
fixing eccentricity properly. The lack of measured eccentricities 
on our other targets shows how difficult it is to measure a small 
eccentricity for these planets as long as no secondary transit is 
detected to constrain it. Spurious eccentricities tend to appear 
in fits to data sets where the radial velocities are not sampled 
uniformly around the orbit, and where the amplitude is small 
compared to the stellar and instrumental noise levels. 

A good example is the case of WASP- 17b for which the dou- 
bling of high precision RV points solely permitted us to place a 
tighter constraint compared to |Anderson et al.| ( |2010| l. 



5.2. Fitting the Rossiter-McLaughlin effect 

Our observations yielded results from which five sky-projected 
spin-orbit angles /3 have been determined with precision better 
than 15°. Three of these angles appear to be retrograde: half our 
sample. Adding the two other stars from our original sample that 
have been published separately (WASP-6b and WASP-8b) we 
obtain 4 out of 8 angles being not just misaligned but also over 
90°. The precision on the angle depends mostly on the spec- 
troscopy as is shown by comparing fits where parameters con- 
trolled by the photometry are kept fixed (in the appendices). 

The error bar on WASP-4b's f3 is large. A degeneracy appears 
when the impact parameter is close to between V sin / and f3. 
The estimate of the spin-orbit angle therefore relies on a good 
estimate of the stellar rotational velocity as well as with getting 
a stronger constraint on the shape of the Rossiter-McLaughlin 
effect. WASP-5b, WASP- 17b and WASP- 18b are also affected 
by this degeneracy, with much lower consequences, when the 
MCMC takes a random step in low impact parameters. This is 
controlled by the use of a prior on V sin /. 

When the planet is large compared to the parent star, or 
the star rotates rapidly, the cross-correlation function develops 
a significant asymmetry during transit. This happens because 
the spectral signature of the light blocked by the planet is par- 
tially resolved. Fitting a Gaussian to such a profile yields a ve- 
locity estimate that differs systematically from the velocity of 
the true lig ht centroid.|Winn et al.| ( f2005| l, and later Triaud et al. 
(2009} and |Hirano et al.| ( |2010|l showed how this effect can lead 
to over-estimation of V sin /. Hirano et al.| ( |2010[ ) have devel- 
oped an analytic method to compensate for this bias. Cameron 
|et al. (2010 1 circumvent the problem altogether by modelling 
the CCF directly, decomposing the profile into a stellar rotation 
profile and a model of the light blocked by the planet. 

Only one star in our sample suffers from this misfit: WASP- 
1 8b where easily we see that the value the fit issues for the V sin / 
is above the estimated value taken via spectral analysis. WASP- 
17b is the second fastest rotating star. If affected, it is not by 
much: the fitted V sin / is found within 1 <x of the v sin / . 



As shown in Fabrycky & Winn (2009 ), we can get an idea of 
the real angle if/ from f3 by using the following equation, coming 
only from the geometry of the system: 



cos if/ = cos / cos i + sin / sin i cos fi 



(2) 



where / is the inclination of the stellar spin axis and i the incli- 
nation of the planet's orbital axis to the line of sight. 

Using the reasonable assumption that the stellar spin axis 
angle / is distributed isotropically, we computed the above equa- 
tion using a simple Monte-Carlo simulation to draw a random 
uniform distribution in cos /. We also inserted the error bars on i 
and /3, using a Gaussian random number adjusted to the 1 <x error 
bars printed in table [5] Fig. [8] shows the transformation from fi 
to iff for our targets, also illustrating the importance of including 
error bars in the calculation. We computed the lower ifi (at the 
3 a limit) and found that in the stars we surveyed: WASP- 15b 
is > 90.3° and WASP- 17b > 91.7° therefore retrograde, while 
WASP-2b is > 89.8° most probably retrograde. 

Statistically we will fail to detect a Rossiter-McLaughlin ef- 
fect (hence f3 and if/) on stars nearly pole-on (with a low /). 
WASP-2b, with its small V sin I could be a close case. It could 
be one reason why its RM amplitude is so small (or stellar rota- 
tion so low). We observe that the spread in if/ is larger than for 
our other targets. 
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0.02 




Fig. 8. top Smoothed histogram of the if/ distribution for WASP- 
5b. The dotted line is when errors on i and (3 are set to zero. 
The plain curve shows the same conversion from (3 to \fi but with 
all errors accounted for. bottom: 6 smoothed histograms of the 
distribution in if/ our six targets: a)-WASP-2b b)-WASP-4b c> 
WASP-5b d)-WASP-15b e)-WASP-17b f)-WASP-18b. Bins are 
of 1°. 



X reduced ^ ^ * 

It can be remarked from the text or from the appendices that a 
few of our objects have Adduced > m me case °f WASP-5 
notably. This shows the models are not adjusted perfectly to the 
data. As showcased by model fits to WASP-4, 15 and 17 (with 
^reduced < an d m many publications using the CORALIE and 
HARPS spectrographs, produced error bars on individual radial 
velocity points are well estimated and understood and worth us- 
ing as they are ( Lovis et al.|2006} . 

An easy way to solve the problem would be to scale error 
bars so as to achieve an acceptable^ 2 . But increasing error bars 
is not making the model more just; rather, it hides that we do 
not understand everything: that an extra signal is observed by 
the instruments; there is information in a bad^ 2 . Error bars can 
be scaled with a value of stellar jitter added quadratically to in- 
dividual errors, but this applies only if one samples randomly 
over long periods of time an extra stellar signal. In our case, part 
of our out-of-transit RVs would feel this jitter, but it would not 
apply in the same way to the Rossiter-McLaughlin sequences 
during which we are sensitive to a correlated noise of a different 
frequency. This renders the increase of error bars prone to errors 
of judgement, thus leading to a wrong computation of the model. 



Hence, we decided to produce results without interfering 
with the way the data is estimated giving the best optimisation 
of the data that we could produce using known and substantiated 
physics. We leave to the reader the assessment of where this ex- 
tra signal originates from. 

5.4. Correlations between parameters 

We present a compilation of results from all known observations 
of the Rossiter-McLaughlin effect in transiting exoplanetary sys- 
tems in Table [5] No clear correlation is evident between impor- 
tant planetary parameters such as radii, masses, eccentricities, 
orbital periods, (5 and V sin /, except that planets with M < 2 Mj 
and e > 0. 1 are rare among transiting systems (the only two are 
Neptunes around M dwarfs); this remark is independent from 
having a Rossiter-McLaughlin measurement or not. It is hard to 
see if this is really a result, or a bias due to observations (eg: 
transits harder to extract from the survey photometry, or to con- 
firm via radial velocity), or a lack of precision during follow-up 
making eccentricity hard to detect with confidence. WASP- 17b, 
for example, was previously thought to be the most eccentric 
transiting planet with M < 2 Mj but our analysis yields only an 
upper limit e < 0.110. Eccentricities with as great as e = 0.1 
have been published for some planets with masses less that 2Mj, 
none of these results are significant at more than the ~ 2 cr level. 

The current (M p sin i, e) distribution in radial velocity does 
not show this result, but these masses are only minimum masses. 

Amongst planets where eccentricity is firmly detected, four 
out of seven are misaligned. Some of the hot Jupiters appear to 
be in multiple systems but this appears unrelated to other pa- 
rameters such as eccentricity or misalignment. Examples are: 
HD 80606 (N aefet al.||2001) , HD 18 9733 (|Bakos et al.|[2006] >, 
Hat-P-1 (Bakos et al . | 2007| ), WASP-8 ( [Queloz et al.|2010|>, Hat- 
P-7 (|Winn etal.|2009b| ), WASP-2, TrES-2 and TrES 4 ( |Daemg"enl 
|et al.|2009| >. ~ 



6. Discussion 

After a long sequence of closely-aligned planets ( Fabry cky &| 



Winn 2009|l, the sudden appearance of so many misaligned plan- 



ets is somewhat surprising if not unpredicted. In a collapsing 
gas cloud, conservation of angular momentum will create a disc 
from which a star can form. Thus it is expected that star and 
disc rotate in the same direction with parallel spin axes. If plan- 
ets form in and migrate through the disc, we can extend the idea 
that planets' orbital axes and stellar rotation axes ought to be par- 
allel. Tides alone cannot make a planet retrograde dHut][l981 



Therefore it is expected that the creation of retrograde planets 
involves another body: planetary or stellar. Several papers (|Wu| 



l et al.|2007||Fabrycky & Tremaine|2007| JNagasawa et al.) 2008 
|chatterjee et al.|2008[|Juric & Tremaine|2008[ |Bate et al.|2000 



2009]) produce via various processes, orbits which are not copla 
nar with the host star's equator. Of these papers|Wu et al. ( 2007 1, 
Fabr ycky & Tremaine| ([2007 ) and Nagasa wa et al.| ( |2008) > pro- 
duce the largest range of angles. 

When combining the 26 RM effects that have been observed, 
we now see that eight planets are sever ely misaligned : XO-3b 
(]Hebrard et al.|2008HWinn et al.|2009c|> , HD 806 06b ( |Moutou| 
|et al.||2009||Pont et al.||2009| |Winn et al.||2009a|i WASP- 14b 
Pohnson et al.|2009|l, Hat -P-7b (IWinn et al.|2009b| |Narita et al.| 
|2009| ), WASP-8b ( |Queloz et al.||2010| i and WASP-2b, WASP- 
15b and WASP17b. Of these eight, five have been found to be in 
retrograde orbits, four from our survey. 
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P(°) 



Fig. 9. top: Histogram of all the j3 measured, binned by 20°. bot- 
tom: Cumulative probability function f or models by |Fabrycky| 
& Tremaine| ( |2007| l (blued dashed) and Nagasawa et al. ( 2008 ) 



(red dotted) converted from ip to fl, compared with current obser- 
vations of /? (plain black). The vertical black dotted line shows 
f3 = 30°. Above that, planets are considered misaligned. 



Three additional targets may be misaligned: Kepler-8b 
( [Jenkins et al.|2010|> CoRoT-lb ( |Pont et al.|20T0] > and CoRoT- 



3b ( Triaud et al. 2009 1. All three are around faint stars and fairly 
fast rotators making it hard to determine the angle. All f3 mea- 
surements have been plotted in Fig. [9^. Because we only mea- 
sure the sky-projection of the angle, the planets can in fact be in a 
variety of configurations. What is their real angle if/ distribution? 

Deconvolving the whole (3 distribution into if/ is hypothesis 
dependent. Hence, to compare the observational data and theo- 
retical predictions we chose to produce cumulative histograms of 
observational and theoretical yS angles in Figj9|). We transformed 
predictions from Fabrycky & Tremaine ( 2007| > and Nagasawa 
et al. ( 2008| l, by taking their if/ histograms and transforming them 
geometrically into observable /3, with the assumption that / is 
isotropic. For a fixed if/, we define an azimuthal angle a mea- 
sured from a zero point where the star's north pole is tilted to- 
wards the observer. If we precess the star for a e [0, 2n[ we 
obtain fi via a Monte Carlo simulation from solving: 



tan ft =: tan if/ sin a 



(3) 



using the conservative assumption that i = 90° since these sys- 
tems are transiting. 



Results from this transform are in Fig.|9j5. The observational 
data has been overplotted. Both observations and models by 
Fabr ycky & Tremaine| ([2007 ) agree that about 55 % of planets 
should appear with ft < 30^ Overall the theoretical distribution 
is a little steeper than the observations. We clearly remark that 
predictions by Nagasawa et al. (2008 ) agree in range but not in 
the shape of distribution of observed /3, notably, it lacks enough 
aligned systems. This model is handy to illustrate the difference 
between the observations and a distribution isotropic in if/. 

Disc migration models would only produce a steep distribu- 
tion reaching unity before 30°. A combination of several models 
is not attempted here because of the vast amount of possibilities 
and the likelihood that models will evolve. 



The theoretical if/ distribution published by Fabrycky & 
Tremaine ( |2007| l, transformed into ft, shown along the angle dis- 
tribution obtained from observations, in Fig. [9J5 gives a remark- 
ably close match. If the form of this distribution is borne out by 
future observations, we may then conclude that hot Jupiters are 
formed by this very mechanism. Wu et al. ( 2007) predicti ons are 
essentially the same as those from Fabrycky & Tremaine (2007 ). 



We also attempted to generalise the method explained in sec- 
tion 5.2 to all objects presented in table [5] we are going to as- 
sume two distributions for the stellar spin axis in order to derive 
a distribution of real obliquity iff. Two hypotheses were checked. 
The first was to assume an isotropy of the stellar axis angle / by 
taking a uniform distribution in cos / from to 1 ; the second was 
to assume stellar axes are aligned with the plane of the sky. For 
this last hypothesis we assumed all cos / followed a Gaussian 
distribution centred on with a variance of 0.1, an error bar cor- 
responding to the best of what observations can give us at the 
moment to constrain the stellar /. Taking these hypotheses al- 
lows to test for extremes and get an idea of the true proportion 
of misaligned planetary systems. 

The first hypothesis is shown on Fig. [T0|plotted in compar- 
ison with the theoretical predictions by Fabrycky & Tremaine 
(2007). Inferring an isotropic distribution in stellar axes gives 
us as an upper bound that 86.2 % of the probability density dis- 
tribution is at t^r > 30°. The other hypothesis gives a propor- 
tion of 43.6% of misaligned systems. The effect of constrain- 
ing the stellar / makes every individual contribution narrower in 
range. Taking a stricter constraint does not change this propor- 
tion much. 

Both hypotheses are at the extremes of what the real distri- 
bution of / is. We interpret these results as showing that between 
45 and 85 % of systems are misaligned with if/ > 30°. Aligned 
systems are no longer the norm, radically altering our view on 
how these hot Jupiters formed. 



Fabrycky & Tremaine 


(2007 


I and Wu et al. 


(2007) 


Kozai mechanism dKozai 


1962; 


Wu & Murray 


2003) 



use the 



by an outer binary companion to the inner planet, to move the 
planet from the ice line where it is thought to form, to the in- 
ner stellar system. As the planet gets closer to the primary, tidal 
friction helps to break the Kozai cycles and finalise the planet's 
orbital parameters. Their equations are extracted from work by 
Eggleton & Kiseleva-Eggleton (2001 ). The resulting if/ distribu- 



tion extends from 0° up to 150° away from the primary's rotation 



axis (see Fig. 10 1. In this scenario, the planet can be created in 
a binary star system, or around a single star which acquired a 



5 this criterion of misalignement of ft > 30°, is a limit where, with 
current error bars on /3, one can generally have a significant detection of 
a misalignement. 
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Fig. 10. The above histogram transformed into the real angle if/ 
in solid line and smoothed to bins of 1 ° . Red dotted curves show 
key individual objects in order to illustrate some of the features 
of the overall distribution. The blue dashed histogram is the re- 
production of the theoretical histogram published by Fabrycky 
|& Tremaine1 ( [2007l ) and solely plotted over. a)-HD 189733b b)- 
XO-3b c)-HD 80606b d)-WASP-8b e)-WASP-15b f)-Hat-P-7b. 
The black dotted line shows if/ = 30°. Above that, planets are 
considered misaligned. 



companion through interactions in its cluster of origin ( Pfahl & 
Muterspaugh 2006) .|Fabrycky & Tremaine| ( |2007| ) following on 
a paper by |Ma lmberg et al. (2007J, also predict that in multi- 



planetary systems undergoing Kozai cycles thanks to a nearby 
star, the most massive planet would survive the resulting planet- 
planet scattering. Although Kozai cycles are usually associated 
with high eccentricities, we should not be surprised by the pres- 
ence of so many misaligned planets on circular orbits. As simu- 
lated in the case of HD 80606b in |Fabrycky & Tremaine| ( [2007) , 
the Kozai cycle has ended (is responsible the close proximity 
of the planet to the central star at periastron making precession 
dominated by general relativity rather than by the action of the 
third body). The planet appears now in a process of circularisa- 
tion that will take ~ 0.7 Gyr, while its angle if/ remains almost 
constant. 

Nagasawa et al. (2008) model scattering processes between 
planets creating a pair where one planet is on a close orbit and 
the other around 40 to 100 AU which then drives Kozai cycles 
on the inner planet. They also use tidal friction with the star. 
These authors predict with orbits with a wide range distribution 
of inclinations and eccentricities which does not reproduce our 



observations as closely as [Fabrycky & Tremaine (2007]) do. All 



other authors fall short of the wide range of angles that we detect. 
This, however does not mean that the processes they describe do 
not happen in combination with those talked about here. In ad- 
dition we cannot rule out that each of the current theoretical dis- 
tributions will evolve thanks to greater scrutiny of their starting 
hypotheses leading to new simulations. Typically, tidal interac- 
tions between the star and the planet have been understudied. 
Any change in the way tidal processes are treated will alter the 
rate at which planets would realign the stellar spin axis ( |Winn 
|et al.|2010| . New effects are also likely to be imagined such as 
these Kozai cycles between a misaligned planet and a disc pre- 



sented in Terquem & Ajmia (2010 ). 



If the Kozai effect were found to be the dominant process 
leading to the creation of hot Jupiters, there is no reason why 



longer period planets should not have undergone similar cycles. 
The only difference would be that having greater periastron dis- 
tances, tidal friction was less active. It would then be expected 
that lone Jupiters on large eccentric orbits be misaligned as well. 
HD 80606b would be part of that population. We could then have 
a lone Jupiter population of which hot Jupiters are a subset, and 
another planet population where dynamical interactions and tidal 
migration did not act. 

7. Conclusions 

The observations reported here bring the total number of transit- 
ing planets with known sky-projected obliquities from 20 to 26. 
Among this enlarged sample, eight show significant projected 
spin-orbit misalignments; and of these eight, five show apparent 
retrograde motion. This projected angle (3 can be transformed 
statistically into the real spin-orbit angle if/. Although 1/3 of 
planets have j3 + 0°, the distribution in if/ shows that up to 85% 
of hot Jupiters are misaligned. The angle range and shape of the 
overall iff distribution appears consistent with the predictions of 
models by |Fabrycky & Tremaine| ( |2007| ) and Wu et al. ( 2007) 1 
using the Kozai mechanism to make planets move inwards and 
tidal friction to reduce their semi-major axis and eventually, cir- 
cularise them. 

Our evidence therefore points towards a dynamical - not lim- 
ited to Kozai - and tidal origin for making hot Jupiters so close to 
their host star. This evidence is the strongest yet to suggest that 
processes others than type I or II migration (using exchange of 
angular momentum between a planet and a disc) are responsible 
for the creation of hot Jupiters. Disc migration alone cannot ex- 
plain the observations; we need to invoke another process. Our 
interpretation is supported by other facts such as how different 
hot Jupiters are spread in semi-major axis compared to multi- 
ple systems (Wri ght et al.| [2009), on how lonely hot Jupiters 
are, and the rarity of hot Jupiters at orbital distances less than 
two Hill radii from the star ( |Ford & R asio 2006 ). These results 
and conclusions should also be a call to account for environmen- 
tal effects on planetary systems in planet formation simulations. 
These systems are not in isolation and interact with their neigh- 
bours. 

We are seeing the coming of a new diversity in planetary pa- 
rameters, coming after large diversities in mass, period, eccen- 
tricity and radius. The variety of angles transformed into ifi, 
is an indication of the physical processes that happened before, 
during and after planet formation. Once again the measurement 
of a new observable has brought a large variety of values reflect- 
ing how rich nature is. 

As more transiting systems are discovered in wide-field sur- 
veys, and follow-up observations of the kind reported here are 
made, the statistical picture that is beginning to emerge will be- 
come clearer. 
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Appendix A: Comparative tables for each star 

Here, for transparency, are the tables recording the results from 
the various fits that were done for each star, which, par com- 
paring them, led to the choice of our solutions, x 1 have been 
tabulated only for the radial velocity data that was used for our 
analysis. In addition, to show where the most important contri- 
butions come from, x 2 have also been added for each set of ra- 
dial velocities separately, as they are presented in the Journal 
of Observations in the following appendix. Finally, because our 
aim was to measure /3, a line with the^- 2 only during the Rossiter- 
McLaughlin effect has been added. Comparisons between the^ 2 
contributions of the overall reflex motion of the star with contri- 
butions during the Rossiter-McLaughlin will show that we tend 
to fit better during transit than outside. The number of points 
during the Rossiter-McLaughlin effect have been chosen as all 
the points measured during transit, plus one point immediately 
on either side when available. 

Appendix B: Journal of Observations 

The Radial- Velocity data extracted by fitting a Gaussian function 
on a Cross-Correlation Function resulting from comparing the 
spectra with a mask corresponding to its spectral type. The data 
is presented per instrument and separated in various datasets: 
overall Doppler shift and Rossiter-McLaughlin effect, as was 
done for the fits. Within each dataset, it is presented chronologi- 
cally. The data is available for download at CDS - Strasbourg. 
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Table A.l. Differences between fits of WASP-2b, 4b & 5b. Adduced nas b een estimated for the radial velocities only. 
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Sharps, rm, n rvs 30.8 ± 7.8 82.7 ± 12.9 53.9 ± 10.4 29.0 ± 7.6 79.1 ± 12.6 51.1 ± 10.1 



WASP-4b 



V sin / Prior 


on 


off 


on 


off 


Fixed Photometry 
on 


V sin / (kms -1 ) 
PC) 

e 

con 


7 14+0.38 

^0.35 
_4+43 


4!f 

4+84 
^-80 


9 i r+0.45 
z - 1J -0.39 

o. +34 

olW 00036 

"• U1UJ -0.0072 

-108+2f 


78+* 1 
28 + i 18 

0106+ 00038 

u -" luu -0.0074 

-107! 280 


2 1 9 +0 ' 35 

c+39 
J -38 



all 56 RVs, 2 sets 

xl w 77.8 ± 12.5 78.0 ± 12.5 75.3 ± 12.4 75.3 ± 12.3 77.8 ± 12.5 

N iof 46 46 44 44 46 

deduced 1.69 ± 0.27 1.70 ±0.27 1.71 ±0.28 1.71 ±0.28 1.69 ±0.27 



/^CORALIE,24RVs 

Sharps, 32 RVs 



28.4 ±7.5 29.1 ±7.6 27.6 ± 7.4 28.0 ±7.5 28.6 ±7.6 
49.4 ±9.9 48.9 ±9.9 47.8 ±9.8 47.3 ±9.7 49.2 ± 9.9 



*harps,rm,15RVs 12.6 ± 5.0 12.5 ± 5.0 12.2 ±4.9 12.0 ±4.9 12.7 ± 5.0 
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Triaud et al.: Six Rossiter-McLaughlin effects observed on WASP planets 
Table A.2. Differences between fits of WASP-5b & 15b. Af^d uce d nas b een estimated for the radial velocities only. 



WASP-5b 



V sin / Prior 


on 


off 


on 


off 


Fixed Photometry 
off 


Vsinl (kms -1 ) 
e 


o 74+0.35 

19 i+lO.O 
1Z - X -8.0 

- 


o 94+0.34 
J - -0.35 

-12.4! 1 - 
- 


Q Q9+0.30 
J ' J -0.32 

_14 i+ro.8 

0.0209+!-!l!!?l 
_ 1 37+14 


3 36 +032 

-i6.r 142 
0.0209+!? , , 7 ! 

w v/ w -0.0087 

-137+[2 


o i q+0.26 

-12.0* 
- 


all 49 RVs, 2 sets 












Xrw 


143.7 ± 17.0 


144.3 ± 17.0 


136.8 ± 16.5 


136.7 ± 16.5 


145.1 ± 17.0 


N iot 


39 


39 


37 


37 


39 


X reduced 


3.69 ± 0.43 


3.70 ± 0.44 


3.70 ±0.45 


3.70 ± 0.45 


3.72 ± 0.44 


/^CORALIE. 16 RVs 

/Sharps, 33 RVs 


20.4 ± 6.4 
123.3 ± 15.7 


20.5 ± 6.4 
123.8 ± 15.7 


26.3 ± 7.2 
1 10.6 ± 14.9 


26.0 ±7.2 
110.7 ± 14.9 


20.4 ± 6.4 
124.8 ± 15.8 


/Sharps, rm, isrvs 


42.8 ± 9.2 


43.9 ± 9.4 


40.6 ± 9.0 


40.7 ± 9.0 


43.9 ±9.4 


WASP- 15b 


V sin I Prior 


on 


off 


on 


off 


Fixed Photometry 
off 


Vsin/ (kms -1 ) 

yS(°) 
e 

a>(°) 


4 26 +0 27 

^'^"-0.32 

139.8!^ 


4 97+0.26 
H -^'-0.36 

139.6!" 


4 o 7 +0.29 
-0.32 

142.6™ 
043 +M20 


4 36 +0 ' 27 
142.7+ 3 - 3 
043+^ : ^ 22 


4 26 +0 ' 28 
139.7!* 


all 95 RVs, 2 sets 












Xrv 


133.1 ± 16.3 


133.3 ±16.3 


130.3 ± 16.1 


130.1 ± 16.1 


133.1 ± 16.3 


N dof 


85 


85 


83 


83 


85 


y 2 

A reduced 


1.57 ±0.19 


1.57 ±0.19 


1.57 ±0.19 


1.57 ±0.19 


1.57 ±0.19 


/^CORALIE,44RVs 

Sharps, 51 RVs 


53.7 ± 10.4 
79.4 ± 12.6 


53.4 ± 10.3 
79.8 ± 12.6 


53.5 ± 10.3 
76.7 ± 12.4 


54.4 ± 10.4 
75.8 ± 12.3 


53.9 ± 10.4 
79.2 ± 12.6 


Sharps, rm, 33 RVs 


47.3 ± 9.7 


47.3 ±9.7 


46.9 ± 9.7 


46.5 ± 9.6 


47.1 ±9.7 



Triaud et al.: Six Rossiter-McLaughlin effects observed on WASP planets 
Table A.3. Differences between fits of WASP- 17b & 18b. deduced nas been estimated for the radial velocities only. 



WASP- 17b 













Fixed Photometry 


V sin / Prior 


on 


off 


on 


off 


off 


Vsin/ (kms -1 ) 

PC) 
e 


Q QO+0.40 

148.5!*' 

~ 


10 14°' 58 

-0.79 

147.3+j!? 


q qs+ - 45 
y - yJ -0A3 

150.9+" 

O A +34 


10 27 +0 6S 
150.5+ 6 , 1 

r\ a/ s 4-'n rnn 
066 

A C + 30 

45 -77 


10 03+ 063 
147 5+4.2 

-4.0 

- 


nil HA PT/c A cvtv 
till 1 Z.H l\v S, H SclS 














190.1 ± 19.5 


190.4 ± 19.5 


187.3 ± 19.4 


186.9 ± 19.3 


191.6 ± 19.6 


N, , 

iv dof 


112 


112 


110 


110 


112 


™ reduced 


1.70 ±0.17 


1.70 ±0.17 


1.70 ± 0.17 


1.70 ±0.17 


1.71 ± 0.17 


/^CORALIE,49RVs 


47.6 ±9.8 


47.4 ± 9.7 


47.2 ±9.7 


47.7 ±9.8 


47.5 ± 9.7 


% CORALIE, 15 RVs 


15.0 ±5.5 


15.0 ±5.5 


16.2 ±5.7 


16.9 ±5.8 


15.0 ±5.5 


Sharps, 16rv s 


23.6 ±6.9 


23.7 ±6.9 


23.5 ± 6.9 


23.7 ± 6.9 


23.9 ±6.9 


X HARPS, 44 RVs 


103.8 ± 14.4 


104.3 ± 14.4 


100.4 ± 14.2 


98.6 ± 14.0 


105.2 ± 14.5 


^CORALIE, RM, 13 RVs 


9.8 ±4.4 


9.5 ±4.4 


9.9 ±4.4 


10.1 ±4.5 


9.8 ±4.4 


X HARPS, RM, 28 RVs 


59.3 ± 10.9 


59.8 ± 10.9 


59.4 ± 10.9 


58.7 ± 10.8 


60.7 ± 11.0 


WASP- 18b 












Fixed Photometry 


V sin / Prior 


on 


off 


on 


off 


off 


Vsin / (kms -1 ) * 
e 

co(°) 


14 04 +0 73 * 
-11 1 +66 


14.66P*. 


14 67 +0 ' 81 * 
-5 +6 2 
0.0084+°™* 
-92.8 + | 2 


i ^ C7+1.01 # 
1 j.j /_ 069 

-4 0+ 50 
-5.0 
n nnss+ 00( » 

U.UOtD „00010 

-92.1^| 


1 ^ SQ+ - 56 * 
iJ,J -0.57 

_4 2 +4 S 
-92.S+ 2 ;? 


all 60 RVs, 2 sets 












2 

^RV 


283.3 ±23.8 


279.3 ± 23.6 


179.7 ± 18.9 


177.8 ± 18.9 


178.4 ± 18.9 


ATdof 


50 


50 


48 


48 


48 


X reduced 


5.67 ± 0.48 


5.58 ±0.47 


3.74 ± 0.39 


3.70 ±0.39 


3.72 ±0.39 


/^CORALIE, 37 RVs 


132.4 ± 16.3 


131.2 ± 16.2 


69.2 ± 11.8 


66.7 ± 11.5 


67.5 ± 10.4 


Sharps, 23 RVs 


150.9 ± 17.4 


148.1 ± 17.2 


110.5 ± 14.9 


111.1 ± 14.9 


1 10.9 ± 14.9 


•Sharps, rm, i2RVs 


113.7 ± 15.1 


110.2 ± 14.8 


98.7 ± 14.0 


98.0 ± 14.0 


98.6 ± 14.0 



Note 4. * these values are not really V sin / but more an amplitude parameter for fitting the Rossiter-McLaughlin effect. Please refer to section 
treating WASP- 18b 
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